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ABSTRACT 

Relations among fermion masses and mixing angles at the scale of grand uni- 
fication are modified at lower energies by renormalization group running induced 
by gauge and Yukawa couplings. In supersymmetric theories, the b quark and r 
lepton Yukawa couplings, as well as the t quark coupling, may cause significant 
running if tan j3, the ratio of Higgs field expectation values, is large. We present 
approximate analytic expressions for the scaling factors for fermion masses and 
CKM matrix elements induced by all three third generation Yukawa couplings. 
We then determine how running caused by the third generation of fermions affects 
the predictions arising from three possible forms for the Yukawa coupling matrices 
at the GUT scale: the Georgi-Jarlskog, Giudice, and Fritzsch textures. 
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1. Introduction 



Thirteen is the number of independent parameters required to describe the 
fermion sector of the standard model: nine masses, three mixing angles, and a 
CP-violating phase. Although such a plethora of arbitrary parameters is usually 
regarded as a weakness, we could instead view the situation as an opportunity to 
reach beyond the standard model. Because the fermion parameters can take on 
arbitrary values in the standard model, any prediction of these parameters can only 
come from beyond the standard model. Conversely, their observed experimental 
values could provide a clue to new physics. 

A step in this direction has been taken with the discovery of various phe- 
nomenological relations among fermion masses and mixing angles, which could be 
viewed as modern-day Balmer formulae. Two types of relations among fermion 
parameters have been explored. The first type links the masses of fermions within 
the same generation to one another. Such relations result naturally from grand 
unified theories (GUTs) when the fermions belong to a common grand unified 
representation and couple to a single Higgs field; group theory then dictates a re- 
lation between their masses [1-3]. The second type of relation connects Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements with ratios of masses of fermions in 
different generations. These relations arise naturally when certain entries of the 
Yukawa matrices vanish, perhaps as a result of discrete symmetries [4, 5]. 

These various relations among Yukawa couplings are presumably a consequence 
of new physics, and therefore hold at the energy scale of the new physics, e.g., the 
grand unification scale. But Yukawa couplings evolve in accord with the renormal- 
ization group (RG) equations, so relations among them that apply at one scale will 
not necessarily hold at another scale. Therefore, before they can be compared with 
low-energy data, GUT-scale relations must be corrected to account for RG running. 

Most of the running of the Yukawa couplings is induced by the gauge couplings. 
For example, group-theoretic relations between quark and lepton masses at the 
GUT scale are greatly modified at low-energy scales because quarks and leptons 
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have different gauge couplings and their masses therefore run differently [1-3]. If 
they are sufficiently large, however, Yukawa couplings themselves induce further 
running [6], and therefore further modifications of relations between masses [7, 
8]. Yukawa couplings also induce running of the CKM matrix elements [9], which 
are invariant under gauge coupling-induced running. The effect of the top quark 
Yukawa coupling on fermion mass and mixing angle relations in super symmetric 
theories was recently investigated in refs. [10-14]. 

In the standard model, the Yukawa couplings of fermions other than the top 
quark are too small to cause significant running. In supersymmetric theories, 
however, the Yukawa couplings of the other third generation fermions, the b quark 
and r lepton, may be comparable to the t quark Yukawa coupling, even though 
their masses are much smaller. This occurs when the expectation value of the 
Higgs field to which b and r are coupled is much less than that of the Higgs field 
to which t is coupled; that is, when tan/3, the ratio of expectation values of the 
two Higgs fields, is large. In this regime, all three third generation fermions may 
cause significant running. This case has also been investigated recently by several 
authors [10, 11, 15-17], who solved the RG equations numerically. 

In this paper, we would like to calculate the effect of RG running induced by 
the entire third generation of fermions analytically. An analytic result would have 
several advantages over a numerical solution. In addition to enhancing intuition 
about the effects of Yukawa coupling-induced running, it would allow one to see 
transparently how changes in the input parameters affect the predictions, without 
having to re-run the numerical routines for each new set of data. It would also 
simplify the error analysis. Although it is not possible to solve the RG equations 
exactly, we introduce an approximation (good to within a few percent) that includes 
the running induced by all three third generation fermions and that allows an 
analytic solution. We then use this approximate solution to determine the effects 
of RG running on several different sets of mass and mixing angle relations, and on 
the predictions that follow from those relations. 
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The RG running of the Yukawa couplings is logically independent of the fermion 
mass and mixing angle relations because the latter arise from new physics at the 
GUT scale whereas the running from the GUT scale to the low-energy scale only 
depends on the particle spectrum below the GUT scale. From the point of view 
of the RG equations, the only role of the fermion relations is to provide boundary 
conditions at the GUT scale. Because of this, we will be able to analyze the 
RG running of the Yukawa couplings independently of any particular set of fermion 
relations. 

One approach is to evolve the matrices of Yukawa couplings down to the low- 
energy scale, and then diagonalize them to find the fermion masses and mixing 
angles. Many degrees of freedom of the Yukawa matrices are not physical, how- 
ever, because of the freedom to perform unitary redefinitions of the fermion bases. 
Therefore, we instead begin with the RG equations for the smaller set of physical 
parameters, viz., the fermion masses, mixing angles, and CP- violating phase. This 
simplifies the computational task by reducing the number of equations, and allows 
us to deal only with physical quantities throughout. 

In sect. 2, we review the one-loop RG equations for fermion masses and CKM 
matrix elements in the minimal supersymmetric standard model. Adopting a non- 
standard parametrization of the CKM matrix, we obtain explicit RG equations for 
the mixing angles and CP-violating phase. In sect. 3, we introduce an approxima- 
tion that allows us to solve the RG equations analytically, including the effects of 
the entire third generation of fermions. We then analyze in sect. 4 the RG effects 
on three different sets of fermion mass and mixing angle relations that might result 
from new physics at the GUT scale. Sect. 5 contains our conclusions. 
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2. Running Masses and Mixing Angles 



In this section, we review the renormalization group running of fermion masses 
and CKM matrix elements. We also include a discussion of phase choices for the 
CKM matrix so that we can obtain explicit RG equations for the mixing angles 
which parametrize it. 

As noted in the introduction, the renormalization group analysis is independent 
of the new physics responsible for relations between fermion masses and mixing 
angles, and can therefore be applied to different sets of relations. We assume, 
however, that the relations result from some grand unified theory, and therefore 
that the three gauge couplings meet at a single scale fx. This can be achieved in the 
context of the minimal supersymmetric standard model, with the supersymmetry 
breaking scale [x susy between 100 GeV and 10 TeV [18]. We will assume that this 
framework describes physics up to the GUT scale fx. 

The one-loop RG equations for the gauge couplings are 

d 

167T 2 — hxgi = kgf, t = \n/x, (2.1) 



dt 



and have solutions 



gf(fx) g 2 8tt 2 \/xJ 
Between fi susy and the grand unified scale fx, the coefficients are given by 

(61,62,63) = (¥,1,-3). (2.3) 

We choose [x susy = 170 GeV for convenience (close to the top quark mass); our 
results will be rather insensitive to the exact value of /x 8VLSy . Using [15] 

4tt "58.5' 4tt "30.1' tor /x susy - 170 UeV, 

(2.4) 
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we obtain 
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JI= 1.2 x 10 16 GeV. 



(2.5) 



4tt 25.0 

Throughout this paper, an overline denotes quantities evaluated at the GUT scale. 

In the minimal supersymmetric standard model, the charge | quarks couple 
to a Higgs field with expectation value (v/y/2) sin/3, where v = 246 GeV and (3 
is arbitrary; the charged leptons and charge — ^ quarks couple to a Higgs field 
with expectation value (v/y/2) cos/3. The fermion masses come from the Yukawa 
couplings 



^Yuk 



where 



V2 



sin 



p) U L YuU R +( ^= cos p) D L Y D D R + 



\V2 



U 



c 



D = 



s 

W 



V2 



cos/3 E L Y E E R +h.c, 



E = 



W 



(2.6) 



(2.7) 



are the fermion fields, and Yjj, Yd, and Ye are arbitrary complex 3x3 matrices. 
(In this paper, we take the neutrinos to be massless.) These Yukawa matrices obey 
the one-loop supersymmetric RG equations [7] 



16n 2 -Y v = 
at 

I6ir 2 -Y E = 
dt 



- cfgf + Tr(3YuYu^ + 3Y V Y^ + Y D Y D * 



cfg 2 + Ti(3Y D Y D ^ + Y e Y e ^) +3Y D Y D 1 + Y v Y^ Y D) (2.8) 



cfg 2 + Tr(3y D Y D t + Y E Y E ^) + 3Y E Y E ^ 



Y E , 



where 



( c i> c 2? C T) — (if> 3, ^) , (cf, C2, cf) — (^, 3, ^y) , (c^, c|, c|) — (||, 3, 0) , 

(2.9) 

between /i susy and /I. 
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Not all the parameters of the Yukawa matrices are physical. Under an arbi- 
trary unitary transformation on the fermion bases, — > LpFp, Fr — > RpF^ 
(where F = U, D, E), the Yukawa matrix undergoes a bi-unitary transformation, 
Yp — > Lp^YpRp, and the charged current becomes off-diagonal, with mixing ma- 
trix Ljj^Lp). We may perform scale-dependent unitary transformations Lp(fi) and 
Rf([i) on the fermion bases so as to diagonalize the Yukawa matrices at each scale. 
Thus 

Yp(fi) = L F \riY F (»)R F (»), F = U,D, E, (2.10) 

where Yp denotes the diagonalized Yukawa matrix, and 

V(fi) = Lu\n)L D (n) (2.11) 

is the corresponding scale-dependent CKM matrix. 

We now derive RG equations for the physically relevant quantities: the Yukawa 
eigenvalues Yp(fj.) and the CKM matrix V(fi) [9]. The transformations on the right- 
handed fields are irrelevant to the CKM matrix, so we begin by writing Yp in terms 
of Lp only 

Y$0x) = Lf\^)Yf^)Yf\^L f ^), F = U,D, E. (2.12) 

Differentiating eq. (2.12), and using eqs. (2.8), we obtain 

^ (y2) = [y* Lu f tu] + _i_ j 6 Tr^)^ + 6 y4 + yygyty* + ygyfg , 

= [yj, Le^Le] + ^ [2 Tr(3Yj + + erj , 

(2.13) 

where = (dLp/dt). The commutator [Yp, Lp^ Lp] has vanishing diagonal ele- 
ments because Yp is diagonal. Thus the RG equations for the Yukawa eigenvalues 



dt 
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y a follow immediately from the diagonal entries of eqs. (2.13). The remaining 
entries of eqs. (2.13) yield equations for the off-diagonal elements of Ljj^ Ljj and 
Ld'Ljj, as long as there are no degeneracies among the quark masses: 

l y\ ~t~ 

{Lu^Lu) a p = —2 ^2 2 _ 2 V °n yi a,/3 = u,c, t, 

71 ~f=d,s,b V P Va 
^ y2 _|_ ^2 

(L D ] L D ) aP = ^2 S 2_ ^7 *7/?> a ^ & a,P = d,s, b. 

(2.14) 

The diagonal elements of Ljj^ Ljj and Ld^Ljj are noi determined by eqs. (2.13). 
It is easy to see why. Equation (2.12) determines Ljj and Lj) only up to right 
multiplication by a diagonal matrix of (scale-dependent) phases. These undeter- 
mined phases contribute arbitrary imaginary functions to the diagonal elements of 
Lu^Lu and L^Ljj. (The off-diagonal elements are unambiguously determined be- 
cause they receive no contribution from the phases.) We can, however, choose the 
phases to make the diagonal elements of Lu^Ljj and Ljj^ Ljj, which are manifestly 
imaginary, vanish: 

(Lu^Lu) aa = (Ld Lo)aa = by an appropriate choice of phases. (2.15) 
The RG equations for the CKM matrix elements (2.11) are then [9] 



ISvr 2 -^ = 16tt 2 (vL^Ln - L v v) 
at V J a 



y=u,c,t s=d,s,b yp t>5 y =u,c,t S=d,s,b y ° y T 

(2.16) 

as long as we choose the phases to guarantee eqs. (2.15). 

We now neglect the contributions to the running caused by the first and second 
generation Yukawa couplings. If we further assume Uu^d'c ^ Dt an d Ud^Vs^ 
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y\, then the CKM matrix RG equations (2.16) reduce to 



1671 " 2 ^ V afi = Vt £ SP V aSV ] StV t p +yl £ -raV a bV%V 7f3 , E a p = < 

8=d,s,b f=u,c,t I _j 

(2.17) 

The CKM matrix elements V a p are not all independent because of the constraint 
of unitarity. 

We prefer to go a step further than previous treatments by deriving RG equa- 
tions for a set of independent quantities parametrizing the unitary CKM matrix: 
the mixing angles and CP- violating phase (but see ref. [19]). To do so, however, 
we must squarely face the issue of phases multiplying the matrix. We adopt the 
following (nonstandard) parametrization of the CKM matrix 










o \ 


/ S1S2C3 + c^e^ 


C1S2C3 - s^e*^ 


S2S3 \ 


/ 




"\ 












S1C2C3 - ciS2e^ 


C1C2C3 + siS2e^ 


C2S3 




e^ s 







I o 





e i& J 


V -S1S3 


-C1S3 


C3 / 


\ 





e i<t>b J 



(2.18) 

where Cj = cos#j and Sj = sin#j, and all the parameters are functions of the scale 
ji. The middle matrix is chosen to have real elements in the third row and column. 
We cannot automatically eliminate the left and right phase matrices by rephasing 
the quark fields; we have already used that freedom to ensure eqs. (2.15), and those 
equations implicitly determine the functions Q (/-O- We can, however, impose the 
boundary condition 4> a (ji) = 0; i.e., the initial values for the five matrix elements 
V a p in the third row or column may be chosen to be real. One can show, using 
VW = 1, that the RG equations (2.17) for this set of five matrix elements close on 
themselves. Since their initial values are real and the coefficients in the equations 
are real, these five matrix elements remain real (i.e., <f) a {fi) = 0) for all fi. In other 
words, the choice of quark phases that guarantees eqs. (2.15) also implies that the 
phases (f> a (fJ>) in the CKM parametrization (2.18) vanish (in the approximation 
that eq. (2.17) is valid). We then obtain the RG equations for the angles 6{ and 



if Va < yp, 
if Va = yp, 
if Va > Vp- 
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the CP- violating phase by substituting eq. (2.18) into eq. (2.17) 



od o d 00 

16-7T —In tan = — y t sin 03, IQtt — In tan 03 = — y t — y b , 

(J.6 (J.6 

167T 2 ^lntan^ 2 = -y 2 sin 2 3 , 16vr 2 -^0 = 0. 
dt dt 



(2.19) 



The RG equations for the Yukawa eigenvalues [7, 9], 







+ 3^ 2 + y 2 cos 2 2 sin 2 3 , 


2 d , 
16tt — \ny c = 




+ 3y 2 + yl sin 2 2 sin 2 3 , 


167r 2 -^-ln^ = 




+ 6y 2 + 2/ 2 cos 2 


16tt 2 — lny d = 




+ y 2 sin 2 0i sin 2 3 + 3y 2 + 2/ 2 , 


16vr 2 ^lnt/ s = 


-4i 


+ y 2 cos 2 6>i sin 2 3 + 3y 2 + y 2 , 


16vr 2 -^ln|/ fe = 


-4i 


+ y\ cos 2 3 + 6yl + y?, 


167r 2 -^-lny e = 
dt 


-49i 


+ 3y 2 + yl 


16vr 2 A ln2/M = 


-c e i9i 


+ Sy 2 b +y 2 , 


16vr 2 — lny r = 


-c-gf 


+ 3y 2 + Ay 2 , 



(2.20) 



are obtained by substituting eq. (2.18) into eqs. (2.13), again neglecting first and 
second generation Yukawa coupling contributions to the running. 

Although we assumed ?/2 < |/ c 2 < y 2 an d y\ <C y 2 <C y\ in deriving eqs. (2.19) and 
(2.20), we did not assume that the mixing angles 9{ were small. Because the third 
generation mixes with the first two, third generation quarks induce some running 
of the mixing angles 9\ and 2 and the ratios of first and second generation quarks, 
Uu/Uc and yd/Us- The amount of running of these quantities, however, is typically 
quite small because of the smallness of 03 ~ 0.05; they change by less than 0.1% 
from the GUT scale to the electroweak scale if y b , yt ^, 1-5. We will therefore be 
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justified in the following in neglecting terms proportional to sin 2 #3 on the r.h.s. of 
the RG equations (2.19) and (2.20). 



3. Approximate Solutions to the RG Equations 



In this section, we find explicit solutions to the renormalization group equations 
for the fermion masses and CKM matrix elements. To do so, we need to make 
several approximations. In deriving the RG equations (2.19) and (2.20) in the last 
section, we neglected the running caused by the first and second generations of 
Yukawa couplings. We now make the further assumption that the mixing angles 
are small. In this approximation, the RG equations for the CKM matrix elements 
simplify to 




for a(3 = ub, cb, td, and ts, 
for af3 = ud, us, cd, cd, and tb, 



(3.1) 



and the Yukawa eigenvalues satisfy 



167T 2 — lny uc 
at 

167T 2 — lny t 
at 

2 d , 

lQn 2 ^-lny b 
at 

16vr 2 — lny eiM 



4gf + 3|/ 2 + y 2 T , 



ct9i + Qyt + vl 



4i + 3y 6 2 + yl 

cfg 2 + y? + Qy 2 + yl 



Ci9i + 3yl 



cfgf + 3y 2 + % 2 . 



(3.2) 



If we define the scaling factors 



AM = exp ^ j£; cf tfV)d(ln//) 



(3.3) 
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and 



B a (jj) = exp 



(3.4) 



then the solutions to eqs. (3.1) and (3.2) are given by 



V nftB, 1 B J 



V a/3 



, for af3 = ub, cb, td, and ts, 
for af3 = ud, us, cd, cd, and tb, 



(3.5) 



and 



VM = y u A u B?, y d (n) = y d A d BlB T , y e {fi) = y e A e B\B Tl 

y c {n) =y c A u Bf, y 8 (fjL) =y s A d BlB T , yM =V ll A e B$B T , (3.6) 

ytfa) = y t A u B?B b , y b {n) = y b A d B t B 6 b B T , y T (pi) = y T A e B b B^, 

where the overline denotes quantities evaluated at the GUT scale. The A a fac- 
tors encapsulate the running induced by the gauge couplings; the B a factors that 
induced by the Yukawa couplings. 

In the approximation of small mixing angles, the four off-diagonal CKM matrix 
elements involving the third generation all run with the same scale factor B^B^ 1 , 
while the remaining five matrix elements do not run at all, as has been observed 
previously [17, 20, 21]. Also in this approximation, the ratios of first and second 
generation Yukawa couplings, y u /yc, yd/ys, and y e /y^, are invariant under running 
induced by third generation Yukawa couplings, as well as under running induced 
by gauge couplings [17, 20, 21]. These results hold to all orders in perturbation 
theory [21]. 

In order to use eqs. (3.5) and (3.6) to scale fermion mass and mixing angle rela- 
tions from the GUT scale to the supersymmetry breaking scale, we must know the 
values of A a (fi susy ) and B a (/i 8VLSy ). The scaling factors due to the gauge couplings 
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Aa(n) obey 



l6ir 2 -\nA a = -c?g 2 , 



(3.7) 



and are easily calculated using eqs. (2.1) and (2.2): 



= n 



i=l 



9i(v>) 



-cf/bi 



i=l 



8tt 2 



(3.8) 



These factors are given at the supersymmetry breaking scale by 

A u (fi susy ) = 3.21, A d (n SVLSy ) = 3.13, A e (fi susy ) = 1.48, for // susy = 170 GeV, 

(3.9) 

using eqs. (2.3) and (2.5). The scaling factors due to the Yukawa couplings B a (fi) 
obey 



167r 2 ^ln^ 
dt 


= vl 


= ft^lBf 


[Bl] , 




16tt 2 — \nB b 
dt 


= vl 


= vU\Bf 


[B 2 t Bl] , 


(3.10) 


16tt 2 — \nB T 
dt 


= vl 


= fAlBf 


[{B b /B T fBl] . 





The B a are equal to 1 at the GUT scale, and decrease monotonically as one lowers 
the scale. The equations (3.10) do not have an analytic solution. We can obtain an 
approximate solution by setting the factors in brackets equal to 1. The equations 
then decouple from one another, and have the solutions 



BtM« [l+f t KM] 1/12 



ln/i 

KM = ^ J A 2 MdM. (3.11) 



In \x 



-1/12 
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We numerically integrate K a (ii) to find 

Ku(Hsusy) = 8.65, K d (n susy ) = 8.33, if e (/i susy ) = 3.77, for // susy = 170 

(3.12) 

The terms in the brackets in eq. (3.10) are all less than 1, so omitting them 
tends to increase the running of the factors B a . Consequently, the expressions 
in eq. (3.11) are smaller than the exact values of B a at /i — fi SUS y by about 1 or 
2%. As we will see in sect. 4, this approximation tends to exaggerate the effect of 
the running induced by the Yukawa couplings. 

In the limit y r <C yt, the approximate solutions (3.11) reduce to the exact 
result for the running induced by the top quark alone [12, 13] 

B t {n)= [l + y 2 t KM]~ 1/U , B b i <l i) = BM = \. (3.13) 



4. Running Relations 

We examine in this section several different sets of fermion mass and mixing 
angle relations, and the effect of renormalization group running on these relations. 
We assume that physics at the GUT scale dictates certain forms, or textures, for 
the matrices of Yukawa couplings. Different physics leads to different textures. In 
this section, we focus on three different textures: the Georgi-Jarlskog texture, the 
Giudice texture, and the Fritzsch texture. We will not be concerned so much with 
the physics behind these textures, but simply take them as given, and examine the 
relations among fermion masses and mixing angles to which they give rise. 

The relations derived from these various texures hold at the GUT scale, and 
need to be scaled down to low energy to yield predictions for measured parameters. 
We use the spectrum of the minimal supersymmetric standard model to run the 
relations from the GUT scale down to the scale at which supersymmetry is broken, 
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A^susy = 170 GeV. Below /i S usy> the CKM matrix elements do not evolve much, 
but the Yukawa eigenvalues continue to run due to QED and QCD effects. This 
additional running is incorporated in the factors i] a , defined by 

r,a = . (4.1) 

In this paper, m a denotes not the physical mass but rather the running mass of 
the fermion, defined by 

v f sin (3 for a = u, c, and t, 
m a = y a (m a )—l (4.2) 
V 2 L cosp tor a = a, s, b, e, /i, and r. 

The physical (pole) mass of the top quark is then related to its running mass by 

rat- (4.3) 



phys 



1 + ^a 3 (m t ) + 0{al) 
on 



In eqs. (4.1) and (4.2), y a (m a ) should be replaced by y a (l GeV) for the three 
lightest quarks, a = u, d, and s. 

When specific numerical values are required in the following, we will use [22] 



m T = 1.7841+°;°^ GeV, m c = 1.27 ± 0.05 GeV, m b = 4.25 ± 0.10 GeV, 

(4.4) 

for the masses, and [13] 



77„ = 2.17, 77, = 2.16, 77c =1.89, % = 1.47, r/ r = 1.02, (4.5) 

for the QCD/QED scaling factors, corresponding to as(Mz) = 0.111. We will 
generally assume that mt is close enough to /i S usy = 170 GeV that running between 



* The t, b, and r Yukawa couplings continue to induce running down to the scale of their 
masses, but the amount of running is much less than that from Jl to ii susy . 
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the two scales is small, 

T) t « 1. 



(4.6) 



There is considerable uncertainty in the values of the scaling factors (4.5) due to 
the uncertainty in a^(Mz) [15, 16]. Since our results are analytic, it is easy to 
determine the effects of choosing other values for m a and r\ a . 



4.1. The Georgi-Jarlskog Texture 



The first texture for the Yukawa matrices we consider is 





(° 


c 


°\ 




( ° 


Fe t( P 





Yu = 


C 





B 


Y D = 




E 







\o 


B 


A) 




V o 





D 



/ F \ 

Y E = F -3E 

V D ) 

(4.7) 

assumed to hold at the grand unification scale. Georgi and Jarlskog [2] originally 
posited this form for the Yukawa matrices in an SU(5) grand unified theory, and 
Harvey, Ramond, and Reiss [3] used it in an SO(10) model. Recently, a number of 
authors [10, 12, 13, 15, 16] have re-examined this texture in a supersymmetric con- 
text. The relations between the Yjj and Yg matrix elements follow if the charged 
leptons and charge — | quarks belong to the same grand unified representation. 
Entries of the two matrices that are equal in magnitude result from Yukawa cou- 
plings to a Higgs field in the 5 of SU(5) or the 10 of SO(10). Entries differing by a 
factor of —3 result from Yukawa couplings to a Higgs field in the 45 of SU(5) or the 
126 of SO(10). The zero entries of Yjj, Yd, and Ye are due to discrete symmetries 
[4, 5] at the grand unified scale. 

The Georgi-Jarlskog texture (4.7) leads to six relations among fermion masses 
and mixing angles. The eigenvalues of the Yukawa matrices (4.7) obey the SU(5) 
relation [1] 



Vb = j 



(4.8) 
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and the Georgi-Jarlskog relations [2] 



Vs-Vd 



3, 



VdVs 



1. 



(4.9) 



The latter two equations can be combined into 



(Vd/Vs 



9 (Ve/Vy) 

[i - (y d /y s )] 2 [l - (yjy,)] 



2 ' 



(4.10) 



The quark Yukawa matrices Yjj and Yjj are diagonalized by bi-unitary transfor- 
mations Yp = Lp^YpRp, with 



(\ \ 



L U = 



C 3 -S 3 
\0 S 3 C 3 / 



/ c 2 -s 2 0' 



/e^ 0" 



S2 c 2 | , L D = 
\0 1 



/ci -si 0\ 



10 1 si ci 

V o o 1/ \ o o i/ 

(4.11) 

and with Rjj and Ftp equal to and Ljj respectively, modulo diagonal matrices of 
signs to make the eigenvalues y a positive. In the unitary matrices (4.11), Cj = cos 9{ 
and Si = sin^i, with 



tan 2£ 1 = ^ 

y s 



t ail 2 2 = ^, 



y c -y u _ yc 
Vt " Vt 



tan 2 fl 3 = ^-^L ~ gg. ( 4 .i 2 ) 



The unitary transformations (4.11) result in a CKM matrix of exactly the form 
(2.18) 



( sis 2 c 3 + cic 2 e^ cis 2 c 3 - sic 2 e*^ s 2 s 3 \ 
sic 2 c 3 - cis 2 e^ cic 2 c 3 + sis 2 e^ c 2 s 3 



V -S1S3 



(4.13) 



-C1S3 



C3 / 



Therefore, in the approximation that the mixing angles are small, the CKM matrix 
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elements satisfy the Harvey-Ramond-Reiss (HRR) relation [3] 




(4.14) 



which depends on the top quark Yukawa coupling, as well as the relations 



U 7 




Vd _ -i<f> \Vu_ 




(4.15) 



which depend on only the first and second generation Yukawa couplings. The six 
relations (4.8)-(4.10), (4.14), and (4.15) hold at the GUT scale. 

The Yukawa matrices do not retain the form (4.7) below the GUT scale because 
of RG running induced by large Yukawa couplings. Rather than follow the evolu- 
tion of the Yukawa matrices, however, we will determine the effect of the running 
on the relations (4.8)-(4.10), (4.14), and (4.15). We use eqs. (3.5) and (3.6) to 
scale the relations from the GUT scale fx to the supersymmetry breaking scale 
A*susy (Whenever the scaling factors A a and B a are written without an explicit 
scale Li throughout this section, fx = fx susy is understood.) The further running of 
the Yukawa couplings from fx susy to the scale of the fermion masses is included in 
the factors rf a defined in eq. (4.1). Thus, the SU(5) relation (4.8) leads to 



m b _ A d rf h B t B\ 
m T A e 7] T B* 

and the Georgi-Jarlskog relations (4.9) and (4.10) imply 

m e m M _ Alrf e rf^ (m d /m s ) 9 {rff JL m e /'n e m li ) 



(4.16) 



m d m s Afros' [1 - {m d /m s )f [l-(^m e /^)f 

since rf d — rf s . The HRR relation (4.14) implies 



(4.17) 




(4.18) 
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and the CKM relations (4.15) lead to 

\V ub \ jr] c m u 



\V U 



m s V 7] u m c 



\V cb \ V t] u m c 



(4.19) 



The running induced by the third generation of fermions is contained in the factors 
B t , B b , and B T . 

Dimopoulos, Hall, and Raby [12] found that four of the six relations implied 
by the Georgi-Jarlskog texture, viz., the two Georgi-Jarlskog relations (4.9) and 
the two CKM matrix element relations (4.15), are unaffected by the running in- 
duced by the top quark, in the limit that the mixing angles are small. We see from 
eqs. (4.17) and (4.19) that, not surprisingly, these four relations are also insensitive 
to the running induced by the bottom quark and r lepton. Therefore, the predic- 
tions of ref. [12] for m s , rrid, <p, and \V u b/V c t,\ remain unchanged even when all three 
third generation Yukawa couplings contribute significantly to the running. 

The SU(5) and HRR relations, however, are modified by the running induced by 
the third generation of fermions. The scaling factors B t) B b , and B T are determined 
by the RG equations (3.10), which have the approximate solutions 

S t w [1 + f t K u ] , B b » [1 + f b K d ] ~ 1/12 , B r » [1 + f T K e ] " 1/12 , 

(4.20) 

where K a = K a (fi susy ) are given in eq. (3.12). The scaling factors B a depend on 
the GUT scale Yukawa couplings y t) y h) and y T . These couplings are related to the 
fermion masses by 

mT = ^^y T BfB 4 T (cosP), (4.21) 
m h = V -^y b B t BlB T { C os{3l (4.22) 
mt = y t B?B b (sinP) « ^B by Jl - B? (sin/5), (4.23) 

using eqs. (3.6), (4.1), and (4.2), where the last equality in eq. (4.23) depends on 
the approximation (4.20). 
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The three GUT scale Yukawa couplings y t , y b , and y T are not independent, 
however. First, the Georgi-Jarlskog texture dictates that y T = y b . Second, from 
the SU(5) relation (4.16), we have 

(The deviation of k from unity shows that significant running must be induced by 
the Yukawa couplings for the SU(5) relation to be valid.) Equation (4.24), together 
with eq. (4.20), may be used to determine y t in terms of y b . Hence, B t , B b , and 
B T , as well as the fermion masses, may be written in terms of a single GUT scale 
parameter y b . This implies a relation between f3 and m t . The parameter (3 may 
be expressed in terms of y b as 

sec/3 = V -^{y h BlB% (4.25) 

using eq. (4.21) and y T = y b . The top quark mass is given by 



m t » ^=B b ^l - fc 12 (B T /B b ) 36 (sin/3), (4.26) 
V2A M 

using eqs. (4.23) and (4.24). Plotting (4.25) and (4.26) parametrically as functions 
of y b , we obtain the relation between mt and tan (3 shown by the solid line in fig. 1. 
To see more explicitly the dependence of m t on f3, we neglect terms of 0(y b ) on 
the r.h.s. of eq. (4.25) to obtain 

sec 3 

m - t£. (4.27) 

then expand eq. (4.26) in terms of y b to obtain the approximate relation 

m t w (185 GeV)(sin/3)[l - (5 x 10" 5 ) sec 2 8 + . . .}. (4.28) 

The Georgi-Jarlskog texture implies an upper bound on the top quark mass, mt < 
185 GeV, which is saturated for tan 8 ~ 10. 
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If the running induced by the b quark and r lepton is neglected, mt increases 
monotonically with /3, as shown by the dotted line in fig. 1, and eq. (4.28) reduces 
to the linear relation between mt and sin/3 found in refs. [12, 13]. The inclusion 
of the effects of all three third generation Yukawa couplings is responsible for the 
last factor in eq. (4.28), and causes m t to decrease for large values of tan (3. Thus, 
a given value of mt may correspond to two possible values of (3 [15, 16]. 

The relation between mt and tan/3 obtained by numerically integrating the 
RG equations (3.10) without approximation is shown by the dashed line in fig. 1, 
and is similar to plots shown in refs. [10, 11, 15, 16]. The solid line, based on 
the approximation (4.20), has the same qualitative behavior as the dashed line, 
but the effects of the Yukawa- induced running for large tan (3 are exaggerated, as 
anticipated in sect. 3. 

The HRR relation (4.18) determines the dependence of |V^,| on j3. Using 
eq. (4.18) together with eq. (4.24), we find 

^ I A e rj T rj t m b m c 
V A d r] b r] c m T m t 

(When running induced by the b quark and r lepton are neglected, this reduces to 
eq. (22) of ref. [12].) Using eqs. (4.25) and (4.29), we plot \V c b\ as a function of 
tan (3 in fig. 2 (solid line). As before, the dashed line indicates the result based on 
numerical integration of the RG equations (similar to plots in ref. [16]). Expanding 
eq. (4.29) in terms of y b , we obtain the approximate relation 

\V cb \ » -^S|[l + (7 x 10- 5 ) sec 2 P + ...]. (4.30) 
Vsm (3 

The Georgi-Jarlskog texture implies that \V cb \ must be > 0.053. 

In principle, given \V c t\, we could use eq. (4.29) to determine (3 (up to a two- 
fold ambiguity) and therefore mt. The uncertainty in |U C &|, however, allows us only 




(4.29) 
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to set bounds on /5. For example, requiring (V^l < 0.058 leads to a lower bound 
on tan/? [12, 13] 



\V cb \< 0.058 =>• tan/5 > 1.5 =>• m t > 155 GeV, 



(4.31) 



as well as an upper bound [15, 16] 



\V cb \< 0.058 tan/5 < 40 =>- m t > 170 GeV. 



(4.32) 



It is the running induced by the b quark and r lepton that causes \V cb \ to increase 
for large tan/5 (4.30), and therefore allows us to derive this upper bound on tan/5; 
neglecting this running leads to the relation shown by the dotted line in fig. 2, 
which obeys \V cb \ < 0.058 for all tan/5 > 1.5. 

Finally, we plot |V^,| as a function of mt in fig. 3. Here our analytic approx- 
imation (solid line) almost exactly coincides with the numerical solution (dashed 
line). Each value of jV^I > 0.053 corresponds to two different values of m t [16]. 
The lower branch of the curve is approximately described by the relation 



and holds for tan/5 < 8. The upper branch of the curve holds for tan/3 > 8; here 
the relation between \V cb \ and m t is approximately linear 




(4.33) 




(4.34) 



which we obtain by combining eqs. (4.28) and (4.30). 
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4.2. The Giudice Texture 



Next we consider the Giudice texture [14] for the Yukawa matrices, 
/0 B 



Y V = 



B 
\B A , 














p e -i<t> 


E 


D 




V o 


D 


C 



/ F \ 
F -3F £> 



\ o d c/ 

(4.35) 

at the GUT scale. (Giudice additionally imposes the ad hoc relation D = 2E, 
but we leave these matrix elements unrelated.) This texture leads to six relations 
among fermion mases and mixing angles. The eigenvalues of the Yukawa matrices 
obey not only the Georgi-Jarlskog and SU(5) relations 



Ve 

Vd 



1 

3' 



^A 4 ~ 9 

— — ~ o, 

y s 



but also the geometric mean relation 



Vb _ , 
_ *-■> 

Vr 



(4.36) 



VuVt = Vc- 



(4.37) 



The quark Yukawa matrices are diagonalized by 




with 



j D 



tan 2 ^i = ^, tan 2 

Vs 

The unitary transformations (4.38) lead to a CKM matrix of the form 






-) 


/ci 


-si 


o\ 


C3 


S3 


si 


Cl 





-S3 


C3/ 


u 





1/ 


= Vu 




h 


D 

<--— ' , 




Vt 






Vb 





( -S1S2S3 + CiC2e*^ -C1S2S3 - Si^e*^ S2C3 

S1C3 C1C3 S3 

V -S1C2S3 - Ci^e*^ -C1C2S3 + Si^e^ C2C3 . 



(4.38) 



(4.39) 



(4.40) 
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Thus, for small mixing angles, the CKM matrix elements obey the relations 



\V ub \^J^. (4.41) 




— 5 



y s V vt 



The matrix element |Vc&| = sin 83 remains arbitrary unless additional constraints 
[14] are imposed on the parameters of the texture (4.35). 

The six GUT scale relations (4.36), (4.37) and (4.41) are modified at low en- 
ergies by RG running. Using eqs. (3.5), (3.6), and (4.1), we obtain the relations 



m e 1 A e n e m u A e r) u ■ . [rru 

m^3^£ ? \Vug\fnJ—, (4.42) 



which are not affected by the running induced by third generation fermions, and 
the relations 

— = (4-43) 
2 

m u = — = B t B b , (4.44) 

Vc m t 



V ^w^t v b 

which are. These relations reduce to those given in ref. [14] if we set B b = B T = 
Vt = 1. 

The relation between m t and tan (3 shown in fig. 1 depends only on the SU(5) 
relation (4.43), and therefore holds for the Giudice texture as well as for the Georgi- 
Jarlskog texture. For the same reason, eqs. (4.24)-(4.28) also continue to hold. 

The up quark mass is determined by (3 in the Giudice texture. The geometric 
mean relation (4.44) together with eq. (4.24) yields 

which we plot as a function of tan/5 in fig. 4. Expanding eq. (4.46) in terms of y b , 
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we obtain 

m u = 2,5 ^ [1 + (1.7 x 10- 4 ) sec 2 /? + ...], (4.47) 
sin p 

implying a lower bound on the up quark mass, m u > 2.5 MeV. We plot the relation 
between m u and mt in fig. 5. Using eqs. (4.28) and (4.47), we find that 

m n 185 GeV lA ao . 

(4.48) 



2.5 MeV m t 



for tan /? < 8, and 



for tan/? > 8. The central value of the up quark mass given in ref. [22], m u = 
5.1 ± 1.5 MeV, could correspond either to a small value of tan/? (~ 0.56) with 
m t ~ 90 GeV [14], or to a much larger value of tan/? (~ 50) with m t ~ 130 GeV 
[15]. An up quark in the range 3.6 MeV < m u < 6.6 MeV requires either 0.4 < 
tan/? < 1.0 with 70 GeV < m t < 130 GeV, or 50 < tan/? < 60 with 100 GeV < 
m t < 160 GeV. 

The Giudice texture also determines |V^| as a function of (3. From eqs. (4.24), 
(4.44), and (4.45), we have 

IVU « , (4^50) 

r] c m t Bl 

which is plotted as a function of tan/? in fig. 6. Expanding eq. (4.50) in terms of 
y b , we find 

f) 0099 

\V ub \ « ^V[l + (1.5 x 10- 4 ) sec 2 /? + ...], (4.51) 
sin ID 

implying the lower bound |V^| > 0.0022. The rather weak constraint \ V u i,\ < 0.007 
requires tan/? > 0.3. 
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4.3. The Fritzsch Texture 



Finally, we consider the Fritzsch texture [5, 23] for the quark Yukawa matrices, 



Yu = 



/ C 
COB 
\0 B A, 



( 



Y D = 



Fe 



t<Pi 



\ 



Fe'^ 1 Ee i<l>2 
\ £e"^ 2 D / 



(4.52) 



at the GUT scale. There is no ansatz for the lepton Yukawa matrix. The matrices 
(4.52) are diagonalized by 



(\ 

c' 3 



\ 



-s' 



/ C 2 -S2 0\ 



\0 s' c' / 



s 2 



C2 







1/ 



/e^ 1 
1 
V e" 



\ 



-i(j)2 



no o \ 



77// 77// 



/ci -si 0\ 

si ci 
\o 1/ 



(4.53) 



with 



tan 2 ^ = ^, 



t a n 2 2 = ^, 



tan z 6> 3 



2 77" Vs-Vd ^ Vs 



Vb Vb 



(4.54) 



When the mixing angles are small, the unitary matrices (4.53) lead to a CKM 
matrix whose elements are given by 





W, 



ub 



V 



cb 




Vc 



V 



cb 




(4.55) 



at the GUT scale. There are no relations among the Yukawa couplings. 



2G 



We scale down the relations (4.55) using eqs. (3.5) and (3.6). The relations 
involving only the first and second generation Yukawa couplings 

IK&I . jVcm u 



Wu 



\V cb \ V Vumc 



(4.56) 



m s y r] u m c 

are the same as in the Georgi-Jarlskog texture, and are not affected by the running 
induced by large Yukawa couplings. The relation for |V^,|, however, is given by 



W, 



cb\ 



in the Fritzsch texture. 




r] s m b 




(4.57) 



The relation (4.57) yields a connection between the top quark mass and \V c i\. 
The small experimental value for \V cb \ requires a large amount of cancellation 
between the two terms in eq. (4.57), which seems to imply both 02 ~ and 
a relatively light top quark [23]. If we neglect the running induced by Yukawa 
couplings and require m s > 120 MeV, m b < 4.35 GeV and m c < 1.32 GeV [22], 
we obtain the upper bound for the top quark mass 

\V cb \ < 0.058 m t < 110 GeV, if B t = B b = 1. (4.58) 

The inclusion of Yukawa coupling-induced running implies an even lower upper 
bound for mt for small or moderate values of tan/3, since typically Bt < B b . If 
tan/3 is very large, however, then it is possible that B b < Bt, loosening the bound 
on the top quark, as Babu and Shafi have pointed out [17]. They demonstrated 
this numerically for a specific value of tan j3. We will derive a modified bound for 
mt using the analytic approximation (4.20). In this approximation, the top quark 
mass is given by 



m t w ^^B b ^J\-B\ 2 for tan/3 > 1. (4.59) 

The largest amount of cancellation will occur in eq. (4.57) when B b /B t is min- 
imized. Holding mt fixed in eq. (4.59), we find that B b /Bt is minimized when 
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1\1/12 



0.85, and therefore 



( 



Bb 
B t 



) 



mm 



150 GeV' 



(4.60) 



Substituting this into eq. (4.57), we obtain 



IV^I < 0.058 m t < 140 GeV, 



(4.61) 



a considerably higher upper bound than eq. (4.58). This high value of the top quark 
mass requires a large value of tan/?. From eq. (4.22), and using the approximation 
(4.20), we have 



In the Fritzsch texture, y T is unrelated to y b so B T is undetermined, but 0.8 < 
B T < 1 for < y T < 2. 



Relations among fermion masses and mixing angles at the scale of grand uni- 
fication are modified at lower energies by renormalization group running induced 
by gauge and Yukawa couplings. In supersymmetric theories, the b quark and r 
lepton Yukawa couplings, as well as the t quark coupling, may cause significant 
running if tan f3 is large. 

In this paper, we have analyzed the running of fermion masses and mixing 
angles caused by the entire third generation of Yukawa couplings. We made several 
approximations along the way. First, we assumed a hierarchy of masses between 
generations. In this approximation, we derived explicit RG equations for the mixing 




(4.62) 



so that 



m t = 140GeV =>• sec(3~50B T . 



(4.63) 



5. Conclusions 
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angles (2.19). Next, we assumed that the mixing angles in the RG equations were 
small; the error from this approximation was shown to be less than 0.1%. Finally, 
we made the approximation that the RG equations (3.10) for the scaling factors 
B a decoupled from one another. This allowed us to obtain the analytic expressions 
(3.11) for these scaling factors, which differ from the exact values by no more than 
1 or 2%. 

We then used the approximate analytic expressions for the scaling factors to 
determine how running induced by the third generation of fermions affects the 
predictions arising from the GUT scale Yukawa matrices. We summarize here the 
results for each of the three textures considered in this paper, using the input data 
(4.4) and (4.5). It is easy to determine the effect of choosing different values for 
the input parameters because our results are analytic. 

The Georgi-Jarlskog texture incorporates the SU(5) relation (4.16), which im- 
plies the relation between mt and tan f3 displayed in fig. 1 and given approximately 
by (4.28). The top quark mass is bounded above by ~ 185 GeV, the bound being 
saturated for tan/5 ~ 10. This texture also implies the HRR relation (4.18); this 
leads to the relation between |V^,| and tan/5 shown in fig. 2 and given approxi- 
mately by eq. (4.30). Consequently, \V c b\ must be greater than ~ 0.053. Requiring 
Wdb\ ~ 0.058 implies that 1.5 < tan/5 < 40, and therefore m t > 155 GeV. 

The Giudice texture also incorporates the SU(5) relation and therefore the 
relation between mt and tan/3 shown in fig. 1. This texture also implies the geo- 
metric mean relation (4.44), which leads to the relation between m u displayed in 
fig. 4 and given approximately by eq. (4.47). The up quark mass is bounded below 
by ~ 2.5 MeV. If the up quark lies in the range 3.6 MeV < m u < 6.6 MeV [22], 
then the Giudice texture implies that either 0.4 < tan/5 < 1.0 with 70 GeV < 
m t < 130 GeV, or 50 < tan/5 < 60 with 100 GeV < m t < 160 GeV. Finally the 
Giudice texture implies the relation between \V U }>\ and tan/5 shown in fig. 6 and 
given approximately by eq. (4.51). 

The Fritzsch texture relates jV^j to m t through eq. (4.57). For small or mod- 
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erate values of tan/5, this relation leads to an upper bound mt < 110 GeV, but for 
large tan/3 (~ 50), the top quark could be as heavy as 140 GeV. 

All of these textures seem able to accommodate a top quark mass in the range 
preferred by the analysis of electroweak radiative corrections, at least for some 
value of tan/5. Which one, or whether any of them, accurately describes reality 
remains an open question. There also remains the more fundamental question: 
what underlying mechanism determines the form of the Yukawa matrices at the 
GUT scale? 
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FIGURE CAPTIONS 



1) The relation between mt and tan/5 implied by the SU(5) relation (4.16). 
Solid line: approximate analytic solution. Dashed line: numerical solution. 
Dotted line: b and r induced running neglected. 

2) The relation between |V C &| and tan/3 implied by the Georgi-Jarlskog texture. 
Solid line: approximate analytic solution. Dashed line: numerical solution. 
Dotted line: b and r induced running neglected. 

3) The relation between \V c b\ and m t implied by the Georgi-Jarlskog texture. 
Solid line: approximate analytic solution. Dashed line: numerical solution. 

4) The relation between m u and tan/5 implied by the Giudice texture. Solid 
line: approximate analytic solution. Dashed line: numerical solution. Dotted 
line: b and r induced running neglected. 

5) The relation between m u and mt implied by the Giudice texture. Solid line: 
approximate analytic solution. Dashed line: numerical solution. 

6) The relation between \V u b\ and tan/5 implied by the Giudice texture. Solid 
line: approximate analytic solution. Dashed line: numerical solution. Dotted 
line: b and r induced running neglected. 
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